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Abstract 

We propose coating spherical particles or droplets with anisotropic nano-sized 
objects to allow micron-scale colloids to link or functionalize with a four- 

(N . 

. fold valence, similar to the sp hybridized chemical bonds associated with, 

<N _ 

. e.g., carbon, silicon and germanium. Candidates for such coatings include 

. triblock copolymers, gemini lipids, metallic or semiconducting nanorods and 

(N . 

conventional liquid crystal compounds. We estimate the size of the relevant 



O 
O 

(N 



X 



nematic Frank constants, discuss how to obtain other valences and analyze 



. the thermal distortions of ground state configurations of defects on the sphere. 

o 
o 



Self-assembly or functionalization of micron-sized particles has many potential uses, in- 
cluding particle-based bioassays catalysis @ and photonic band gap materials ||. How- 
ever, attaching a small, predetermined number of chemical groups, DNA strands or kinesin 
molecules to considerably larger colloidal spheres can be an inefficient process, since the 
number of ligands per sphere is highly variable. If the ligand concentration in solution is 
kept very low, many spheres will have no ligands at all. Higher ligand concentrations require 
tedious sorting methods such as gel electrophoresis to select the spheres with the correct 
number of attached objects 0. 

Controlled fabrication of tetravalent colloids would be of particular interest. Dense glassy 
and crystalline colloidal arrays (possibly involving DNA linker elements ||) usually display 
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the large number of near neighbors [Z = 12 — 14) characteristic of an isotropic pair potential 
||. We analyze here a means by which micron-scale colloids could link or be functionalized 
with, say, a four-fold valence, similar to the sp 3 hybridized chemical bonds on an Angstrom 
scale associated with, e.g., carbon, silicon and germanium. We suggest in particular coating 
spherical colloid particles (or polymerizable liquid droplets) with anisotropic objects such 
as metallic or semiconducting nanorods, gemini lipids, triblock copolymers or conventional 
nematogens. If the particles are sufficiently anisotropic, a two-dimensional nematic liquid 
crystal phase will appear on the surface. The spherical topology then forces four strongly 
repulsive disclination defects into the ground state, a circumstance which would allow cre- 
ation of novel tetravalent colloidal materials with chemical linkers or DNA strands anchored 
at the defect cores. In contrast to fee and bec colloid arrays, a tetravalent colloidal crystal 
with a diamond lattice structure and appropriate dielectric contrast is predicted to have 
a very large photonic band gap 0. More generally, colloidal particles with a 1-, 2-, 3- or 
4-fold valence would also allow creation of functionalized micron-sized objects similar to the 
molecules characteristic of organic chemistry, and could be useful for bioassays and catalysis. 

A two-dimensional nematic at a flat interface is described by a Frank elastic free energy, 
namely [|S| 



where n(r) is a slowly varying director field describing the local molecular alignment and 
K\ and K3 are, respectively, splay and bend elastic constants. (The twist elastic constant 
K 2 is absent in two dimensions.) The effect of thermal fluctuations (which are very strong 
in d = 2) in a sufficiently large system is to renormalize these elastic constants to a single 



point-like defects in nematic textures depends sensitively on the symmetry of the constituent 
particles. For tilted molecules (as might be the case in freestanding smectic liquid crystal 



projection of the molecules onto the plane. The 2rc rotational symmetry leads to defects like 
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(1) 



common value K ||. As a first approximation, we can 



take K « J(KiK 2 ). The nature of 



layers [|ItJ or Langmuir-Blodgett films flTT|), the order parameter is a vector representing the 



the vortex shown in Fig. la. An unbinding of vortex-antivortex pairs causes a transition 
out of the ordered state whenever K < 2/^T/tt |12|]. In contrast, the order parameter for 
molecules which lie down completely on a liquid or solid interface has the higher symmetry 
of a headless arrow. The symmetry under rotations by 7r which results would apply to 
nematically ordered monolayer phases of, e.g., gemini lipids ABA triblock copolymers 
14fl and nanorods made of gold |15|], BaCrC>4 flnj or CdSe Hl7| . In this case, it is well known 



18[| that the vortex in Fig. la is unstable to the two disclinations shown in Fig. lb. 



The 180 degree rotational symmetry allows these mutually repelling defects to separate — 
these defects would be bound with a linear potential for an order parameter with a vector 
symmetry. Disclination unbinding leads to a transition at a temperature which is four times 
smaller, i.e., for K < ksT /2n. 
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FIG. 1. (a) Vortex configuration of a vector order parameter in the plane. The order parameter 
rotates by 2tt on a circuit around the defect at the center, (b) The texture in (a) is unstable to 
two disclinations for a configuration of headless vectors in the plane. The order parameter rotates 
by 7r on a circuit around each of the defects marked by small dots. 

Although pairing of defects and antidefects leads to a virtually defect-free equilibrium 
state at high densities or low temperatures at flat interfaces, the situation is quite different 
on surfaces with a nonzero Gaussian curvature. For liquid crystal phases embedded in an 
arbitrary curved surface, the free energy in the one Frank constant approximation reads 

F = \kJ d 2 x^gjx)[d^ + Ti^Y, (2) 

where x = (xi,^) represents a set of internal coordinates, the director is n = n(xi,X2), 
g is the determinant of the metric tensor gij(x), and the {T J ki } are connection coefficients 
associated with parallel transport of the order parameter. For a rigid sphere of radius R 
with polar coordinates (8, 0), we have y/g = R 2 sm8, = — sin 6 1 cos 6 1 , rjL = rt e = cot 8 
and all other Y 3 ki = 0. The free energy (2) also describes hexatic liquids (i.e., liquid crystals 



with a s^r-fold symmetry) on curved surfaces fl9| , |2"0|| . In this case, the Gaussian curvature 
of the sphere produces an unpairable excess of 12 "disclinations" (i.e., defects around which 
the order parameter rotates by 2tt/6 = 60 degrees) at the vertices of an icosahedron in 



the ground state [JTS]] . In an elegant paper, Lubensky and Prost have studied the ground 



states of more general order parameters with a p-fold symmetry on the sphere pT] . In this 



case, the order parameter rotates by 2it/p on a circuit which encloses the minimum energy 
defect, and the j-th defect is characterized by an integer rij = ±1, ±2, • • -, which specifies 
the winding angle 27m.,- /p. The Poincare index theorem then implies that 



5>,- = 2p, (3) 



j 



for any ordered texture embedded in a surface with the topology of sphere. With our con- 
ventions, positive winding numbers are favored by positive Gaussian curvature and negative 
winding numbers are favored by saddle-like regions of negative Gaussian curvature. 

As pointed out by Lubensky and Prost, the ground state of a two-dimensional nematic 
texture on a sphere consists of four rij = +1 disclinations at the vertices of a tetrahedron 
| |2 1| | . Aligned states of vector order parameters on the sphere, on the other hand, will be 
interrupted by just two rij — +1 vortices, corresponding to the familiar north and south pole 
singularities of the lines of latitude and longitude of spherical cartography. Splay dominated 
textures for vector and headless vector order parameters on the sphere are illustrated in 



Figs. 2(a) and 2(b), respectively [23]. The texture in Fig. 2(a) can in fact be transformed 
into those of Fig. 2(b) by invoking the instability shown in Fig. 1. The inevitable defects 
associated with liquid crystalline textures on the sphere could be exploited to create colloids 
with a 2- and 4-fold valence in a number of ways. If, say, gold is deposited on coated solid 
spheres, the texture will act as mask such that the two or four "bald spots" at the defect 
cores will be coated preferentially. One could then attach, e.g., thiol chemical linkers via the 
extra gold at the bald spots. Alternatively, one could introduce DNA linkers attached to 
amphiphillic molecules on polymerizable liquid drops. Linkers immiscible with the aligned 
liquid crystalline molecules will segregate preferentially at the defect cores, similar to impu- 
rity segregation at grain boundaries in conventional metallurgy. While we do not wish to 
minimize the experimental challenges, the examples cited above illustrate the possibilities 
for directional bonding associated with ordered states in spheres. Topology can be used 
to create a micron-scale directional chemical bonds similar to those produced by quantum 
mechanics on an Angstrom scale. 
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FIG. 2. (a) Splay vector order parameter configuration on a sphere. There are +1 vortices at 
the north and south pole, (b) Two views of a splay configuration of headless vectors on the sphere. 
There are four disclinations at the vertices of a tetrahedron. (Note that one row of nematogens 
tracks the seam of a baseball!) 

To explore the fidelity of directional bonds created by the above scheme, we need a more 
detailed analysis of ground states on the sphere and their stability to thermal fluctuations. 
To this end, we exploit an alternative representation of Eq. (2) in terms of a defect density 

ihr^ .„„ ,,r— 



s(x) = —J2 n j s{2 \ x - x j)/y9{xj) 



V 



and the Gaussian curvature G(x) namely |24 

F = \K J ifxyf^j J d 2 y^y)\G{x) - s(x)}^\ [G(y) - S(y)]. 



(5) 
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(6) 



Here, 1/ A denotes the inverse of Laplacian operator, 

vf 3 

appropriate to a curved surface with metric gij(x) and the defects sit at positions {xj} with 
a set of integer winding numbers {%}. The resulting energy for liquid crystal textures with 
a p-fold symmetry on the sphere is quite simple, 



6 



F = 



ttK 



^riiTij ln(l - cos/%) + E(R)^n 2 . 




(7) 



Here the angle fly = dij/R, where dij is the geodesic distance between defects i and j, and 
E{R) is a defect self-energy, 



The quantity E c is a microscopic energy which depends on the details of the particle inter- 
actions near the defect core on length scales of order a, where a is the particle separation. 

The repulsive interaction between defects in (7) is very long-range, extending out to a 
significant fraction of the circumference of the sphere. Unlike conventional Coulomb inter- 
actions between charged particles, this interaction cannot be screened if the order on the 
surface is strong. The defect self-energy (8) diverges logarithmically with the sphere radius 
and accounts for the instability indicated in Fig. I. To see this, note that the 2ir vortex 
of Fig. la in a medium of aligned headless vectors corresponds to a topological charge 
n = 2 with p = 2. Two such defects at the north and south poles contribute a logarithmi- 
cally diverging term in the energy, 2[(n/p) 2 nK\n(R/a)} = 2irK \n(R/a). If, however, each 
defect splits into two disclinations with charge n — 1, the corresponding contribution is 
4[(n/p) 2 nK ha(R/a)] = 'kK\tl{R/ a). We see that the dominant term in the energy for large 
R/a has been reduced by a factor of two. It is easy to check that the ground state energy 
configuration with the repulsive interaction in (7) is indeed two vortices at the north and 
south poles, and four disclinations at the vertices of a tetrahedron for vector and nematic 
ordering, respectively. 

The long-range interaction in Eq. (7) will resist thermal disruption of the ordered defect 
arrays mentioned above. To investigate this effect for vector order parameters, we define a 
bending angle 6 between the two antipodal defects by setting = n + 9. Upon expanding 
(7) in 6, we find that F const. + (tt/4)K6 2 . The equipartition theorem now shows that 
for valence Z = 2 colloids in the limit K » ksT, 



E{R) 



irK 



Y \n(R/a) + E c . 



(8) 



<cosA>~-l + ^ 2 > 
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Determining the thermal disruption of the tetrahedral ground state which results for nematic 
order on a sphere is more involved. To get a rough estimate, imagine freezing three of 
the disclinations (with j = 1 — 3) at the vertices of a tetrahedron and then allow the 
remaining defect (with j = 0) to fluctuate away from its equilibrium position on the z 
axis. Upon expanding Eq. (7) in the polar coordinates (9, 0) of this defect about an 
equilibrium configuration with tetrahedral angles = cos _1 (— 1/3) (z ^ j), we obtain 
F m const. + (3ir/lQ)K9 2 , which leads to (9 2 ) = SksT/SirR. A more precise analysis 
follows from rewriting Eq. (7) (with p = 2 and all rij = 1) as 

F = -— E ln (l-m l -m J )+4E( J R), (10) 

where the nij are three-dimensional unit vectors specifying the positions of the defects on 
the sphere. Upon writing rrt,- = m" + 8m.j, where 8uij _L represents a small deviation 
from a perfect tetrahedral configuration such that m? • mj = —1/3 (i ^ j), the energy to 
quadratic order is 

ttK 

i,j a/3 



F = const . + — E E M^5m?5m^ (11) 



where a and (3 run over the two independent components of the {<5mj} and is an 8 x 8 
matrix describing the deformations of a perfect tetrahedron on the sphere. The eigenvalues 
of this matrix can be classified according to the irreducible representations of the symmetry 
group of the tetrahedron, similar to the analysis of the vibrational modes of the methane 



molecule [25j. A somewhat lengthy calculation shows that the eigenvalue spectrum of M^ 13 



reads 

{X j } = (3/8){0,0,0,l,l,2,2,2}. (12) 

Rigid body rotations generate the three normal modes with zero eigenvalues. The remaining 
eigenvectors generate two- and three-dimensional representations of the tetrahedral symme- 
try group. The doublet corresponds to a shear-like twisting deformation of the tetrahedron 
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while the non-zero triplet of eigenvalues represents stretching or compression of vectors join- 
ing neighboring defects. Upon re-expressing the free energy (11) in terms of the normal 
modes we calculate that the deviation for valence Z = 4 colloids from perfect tetrahedral 
angles due to thermal fluctuations is given by 

1 1 k T 

< cos ^-3 + 97ir {Z = A) - (13) 

Because the distortion from perfect alignment for Z = 2 and Z = 4 colloids is propor- 
tional to ksT/K, it is of some interest to determine this quantity for a nematic phase formed, 
say, from anisotropic rods [15-17]. Close to the nematic ordering transition, we expect the 
universal values [9, 12] Ky — K^ — K — 2k B T j p 2 n , suggesting fairly large fluctuations in the 
bond angles associated with the topology-induced directional bonding for p = 1 {Z = 2) and 
p = 2 {Z = 4). However, these angles are more robust to thermal fluctuations in the high 
density limit where K » k B T . Suppose the rods have length £ and diameter d. Estimates 
are conveniently constructed by adapting results from an Onsager theory for the bulk Frank 
constants k[ 3 ^ and in three dimensions, a theory valid in the limit £ » d \2E\. For a 



(3) 

monolayer film of thickness d, we have Ki pa dK\ . For well developed nematic order, the 
numerical results of Lee and Meyer when applied to thin films of thickness d lead to the 
approximate formula 



K ^ s]K x K z ^ 0.10p 2 tk B T, (14) 

where p is the areal density of nematogens. This density is of order p ~ l/£ 2 at the ordering 
transition but rises to p ~ l/£d deep inside the nematic phase. Thus, in the dense limit 
K = 0.10 (£/d) 2 k B T, so that k B T/K pa 10 (d/£) 2 . We see that the bond angle fluctuations 
become quite small for large £/d, which is of order 10 in the experiments of Ref. 17. 

The Onsager theory calculations of Ref. []26] show that the bend elastic constant K 3 can 
be much larger than the splay constant K\ for highly anisotropic hard rods. This asymmetry 
favors the "lines of longitude" ' texture of Fig. 2a over a texture which follows the lines of 
latitude. These two configurations are degenerate in the one Frank constant approximation. 
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A latitudinal texture would dominate if K3 << K\. The nematic texture which replaces 
Fig. 2b in this limit is shown schematically in Fig. 3, together with possible polymer linking 
elements emerging from defects at the vertices of a tetrahedron. Investigations of these 
textures for arbitrary K1/K3 are currently in progress using a lattice model of a nematic 
Z7| . Direct simulations of two-dimensional hard rod fluids [ 28 1 on a sphere would also be of 



interest, as would investigations of two-dimensional smectic order on a sphere [27]. 




FIG. 3. Bend texture of headless vectors on the sphere appropriate to K\ > K%. The preferred 
configuration for K\ < K3 is the splay texture shown in Fig. (2b). Polymer linkers are shown 
schematically emerging from the tetrahedron of four disclinations. 

We have argued that coating of colloidal particles with vector or nematic degrees of 
freedom could be used to construct micron-sized objects with valence Z = 2 or Z = 4. Other 
valences and mechanisms are also possible. For example, the bulk nematic liquid crystal 
droplets which appear in polymer dispersed liquid crystals typically have two surface defects 
||29|| , which might also lead to particles with functionality Z = 2. Biaxial nematic droplets 
should have a single "boojum" defect on their surface flT8[] , which upon polymerization might 
lead to a Z = 1 micron-sized analogue of the hydrogen atom. Finally, tetravalent colloids 
with, say, DNA linkers could be tethered via one of the four defects to a surface (see Fig. 
4). This configuration could be useful for colloid-mediated catalysis [2], and the remaining 
crosslinks in the plane of the surface would have a 3-fold valence, i.e. a colloidal analogue 
of sp 2 hybridization. 
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FIG. 4. Tetravalent colloids tethered to a planar substrate. The remaining linkers can connect 
with a 3-fold valence Z = 3. 
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FIGURES 

FIG. 1. (a) Vortex configuration of a vector order parameter in the plane. The order 
parameter rotates by 2tt on a circuit around the defect at the center, (b) The texture in 
(a) is unstable to two disclinations for a configuration of headless vectors in the plane. The 
order parameter rotates by 7r on a circuit around each of the defects marked by small dots. 

FIG. 2. (a) Splay vector order parameter configuration on a sphere. There are +1 
vortices at the north and south pole, (b) Two views of a splay configuration of headless 
vectors on the sphere. There are four disclinations at the vertices of a tetrahedron. (Note 
that one row of nematogens tracks the seam of a baseball! 

FIG. 3. Bend texture of headless vectors on the sphere appropriate to K x > K s . The 
preferred configuration for K x < K 3 is the splay texture shown in Fig. (2b). Polymer linkers 
are shown schematically emerging from the tetrahedron of four disclinations. 

FIG. 4. Tetravalent colloids tethered to a planar substrate. The remaining linkers can 
connect with a 3-fold valence Z = 3. 
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